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ABSTRACT 


In  this  paper»  the  two-population  aodel  for  a carrier-borne  epldeaic  posed 
by  Bailey  (The  Mathsaatlcal  Theory  of  Infectious  Diseases  and  Its  Applications. 
1975,  p.  211)  Is  foroulated  In  a aatheaatlcally  tractable  aianner.  This  nodel 


reflects  the  epidemiology  of  diseases  such  as  malaria,  where  the  progress  of  the 
disease  depends  on  the  Interaction  of  a population  of  mosquitoes  and  a population 
of  humana.  An  expresalon  for  the  mean  duration  time  of  the  epidemic  Is  obtained 


and  a computationally  feasible  algorithm  Is  presented.  Results  of  a study  Investi- 
gating the  conaequences  on  the  mean  duration  time  of  varying  the  Infection  and 
removal  rates  In  the  two  populations  are  given. 


KETM0RD8:  carrier-borne  epidemic,  malaria,  mean  duration  time 
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1.  nmtoDocTioii 

la  elM  ehaory  of  opldoale  procMooa  dovalopod  thus  far,  it  haa  gaaarally  baaa 
aaauaad,  la  tha  alaplaat  caaa,  that  tha  populatloa  of  lataraat  la  cloaad  aad  aay 
ba  dlvldad  lato  ai^populatloaa  ahoae  aaabars  aay  ba  elaaalflad  aa  auacaptlbla  or 
laf active.  A allshtly  aora  coaplicatad  aodal  latroducaa  a third  aub-popalatloa 
of  TMovala. 

Thasa  alapla  aodela  aad  the  aaaoclatad  theory  are  laadaquata  la  tha  atudy  of 
aaveral  dlaeaaaa  (aalarla,  for  axanple)  la  which  the  apread  of  the  dlaaaaa  by 
carriara  la  a recognized  phanoaenon.  A carrier  la  an  Individual  who  aay  tranaalt 
tha  dlaaaaa  to  other  Individuals,  but  who  has  no  overt  dlsaasa  syaptoas  hlasalf . 

In  tha  first  significant  work  in  this  area,  Weiss  (1965)  considers  a cloaed 
population  of  a susceptibles,  into  which  n carriers  are  initially  introduced. 

There  ia  no  subsequent  Introduction  of  carriers,  nor  do  any  susceptibles  bacoae 
carriers  during  the  course  of  the  epidemic.  The  carriers  are  detectable  only  by 
the  dlacovery  of  Infected  persons,  and  the  epidemic  progresses  until  either  all 
suscaptiblas  become  Infected  or  all  carriers  are  renaved. 

Several  extensions  to  the  Weiss  model  have  been  considered.  Downton  (1968) 
allows  for  the  further  creation  of  carriers  from  the  susceptibles  during  the 
course  of  the  epidemic.  Diets  and  Downton  (1968)  relax  the  assumption  of  a closed 
population  and  allow  for  immigration  of  susceptibles  and/or  carriera.  Pettigrew 
and  Waiaa  (1967)  consider  an  epidemic  Involving  two  types  of  infectives  (e.g., 
clinically  Infected  and  sub-cllnlcally  infected  indlvlduala)  in  an  infinite  popu- 
lation of  suaceptlbles . Applying  results  from  branching  proceaaes,  they  obtain 
equatlona  for  the  mean  nuari>er  of  clinically  infected  Indlvlduala  and  for  the  mean 
number  of  sub-clinically  infected  indlvlduala  (carriers)  at  any  time  t. 
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ee«e«ni,  liowwwr,  la  not  with  cha  oaa-population  aodal  of  thasa 
auchera.  but  rathar  with  tha  two-population  nodal  aa  dlacuaaad  by  Ballay 
(197S,  p.  211).  Thla  parnlta  conaldaratlon  of  auch  carrlar-borna  apldanlca 
•a  nalarla.  In  which  tha  carrlara  (aoaqultoaa)  forn  a population  dlatlnct 
fron  tha  boat  population  (hunana). 

In  Sactlon  2 of  thla  papar,  wa  praaant  a atochaatlc  nodal  for  the  two- 
populatlon  carrlar-boma  apldenlc.  A Chaoratlcal  axpraaalon  for  the  naan 
duration  tlaa  la  darlvad  In  Sactlon  3 and  a computationally  faaalbla  algorlthn, 
baaad  on  a racuralon  ralatlon,  la  praaantad  In  Sactlon  4.  In  Che  final  aectlon, 
we  praaant  aona  nunerlcal  raaulta  for  aalected  paraneter  valuea. 
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2.  THE  CARRIER-BORNE  EPIDEMIC  MODEL 

Sj(t),  Ij(t)»  and  Rj(t)  denote  the  nuiri>er  of  ausceptlbles,  Infectlvea, 
and  rMovala,  reapectlvely.  In  population  j at  time  t.  Population  J la  of  fixed 

alxe 

Nj  - Sj(t)  + Ij(t)  + Rj(t), 

J ~ for  all  t.  The  aubacrlpt  j ■ 1 will  refer  to  the  boat  population  and 

J ■ 2 will  refer  to  the  carrier  population.  Let  Z(t)  - (S, (t),I, (t) ; S-(t),I.(t)) 

1 1 z z 

have  realization  z(t)  ■ z ■ (aj^,lj^;  a2tl2)*  Then  Z takea  valuea  In  the  atate 
apace 

A - {z  - (Sj^.lj^;  82,12):  0 S Sj.lj  5 N^;  0 5 Sj  + 1^  i N^,  J-1,2}. 

Note  that,  alnce  Rj(f)  " not  Include  the  nuiri>er 

of  removed  Indlvlduala  In  the  definition  of  g(t). 

In  a almpllflcatlon  of  the  epidemiology  of  malaria,  a human  auaceptlble 
acqulrea  Infection  through  being  bitten  by  an  Infective  female  moaqulto.  The 
pntaaltea  releaaed  by  the  moaqulto  multiply  In  the  liver  of  the  human  boat,  and 
then  aettle  In  the  red  blood  cella.  Female  moaqultoea,  aeeklng  the  blood  re- 
<)ulred  for  the  development  of  their  egga,  are  then  Infected  by  biting  an  Infected 
human  boat.  Thua,  the  probability  that  an  additional  human  auaceptlble  becomea 
Infected  during  a abort  time  Interval  la  proportional  to  the  nuiid>er  of  human 
auaceptlblaa,  the  number  of  Infected  carriers,  and  the  length  of  the  tine  Interval. 
Similarly,  the  probability  an  additional  carrier  becomes  Infected  during  a short  time 
Interval  la  proportional  to  the  nuBd>er  of  susceptible  carriers,  the  nuhbcr  of  Infected 
huauins,  and  the  length  of  the  time  Interval.  The  removal  rates  are  asaumed  linear. 
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Thus,  the  probabilities  of  permissible  transitions  In  the  Interval  (t,t-Mt)  are: 

(I)  P{(Sj^,lj^;s2,l2)  (Sj^-l,lj^+l;82fl2)J  “ ^1*1^2***^^  o(dt), 

(II)  P{(8j^,lj;82,l2)  (Sj^,l2-l;82,l2)}  ■ >1*1 

(III)  P{(8j^,lj^;s2.l2)  (8j^,l2;82-ltl2+l))  ■ ^2*2^1***^^  o(dt), 

(Iv)  P{(Sj^,lj^; 82,12)  (sj^,!^; 82,12-1)}  ■ Y2l2<*t  + o(dt). 


(v)  P{(8j^»l2; 82,12)  (82,l2;82,i2)}  “ ^"^®1®1^2  ^ ^1^1  ^ ®2®2^1  ^ ^2^2^**^ 

+ o(dt), 

where  0^  and  the  Infection  and  removal  rates,  respectively,  are  assumed 
constant . 


3.  MEAN  DURATION  TIME 


Our  epidemic  process  Is  said  to  be  terminated  whenever  no  more  human  sus- 
ceptlbles  can  acquire  Infection.  This  occurs  when  one  or  more  of  the  following 
subsets  of  A Is  entered,  viz.. 


the  nuaber  of  human  susceptlbles  Is  zero, 

E2:  the  number  of  Infectlves  In  each  population  Is  zero. 


E^:  there  are  no  susceptible  or  Infective  carriers  (l.e.,  all  carriers 


have  been  removed) . 


Then,  let 


E * E^  u E2  u ^3* 


where,  from  the  definitions  above, 


E^  - {(0,l2;82.l2)}» 


\ 
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Let  the  conplenent  of  the  set  B be 
G - A-B. 

Now,  we  are  Interested  in  T^,  the  mean  time  required  for  an  epidemic  which  Is 
currently  In  state  m e G to  enter  a state  In  the  set  B. 

To  derive  expressions  for  the  mean  duration  time,  we  consider  our  process 
as  a continuous  time  Markov  chain  with  a one-dlmenslonal  state  space.  In  order 
to  Identify  each  vector  In  A by  a natural  ntiiid>er,  we  apply  the  transformation 

k 3 k(.j,l^;,j,ip 


where 


+ (Nj^^+l)Nj^^(lffl)/2  + (WDsj^  + (S2+I). 


Nj^-Sj+lj,J-l,2;  N-max  (N2^.N2^.)  and  6(x)  - I(x  i 0). 


The  transformed  state  space  Is 
Aj  - {h;  1 s h i 


h Integer}. 


For  example,  the  first  nine  states  In  A^,  regardless  of  the  population  sixes, 


s 

k(£) 

(0,0 

0,0) 

1 

(0,0 

0,1) 

2 

(0,0 

1,0) 

3 

(0,1 

0,0) 

4 

(0,1 

0,1) 

5 

(0,1 

1,0) 

6 

(1,0 

0,0) 

7 

(1,0 

0,1) 

8 

(1,0 

1,0) 

9 
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We  note  that  there  Is  a unique  one-to-one  correspondence  between  the 
vectors  In  A and  the  Integers  In  A^.  An  algorithm  to  determine  the  vector 
In  A which  corresponds  to  a given  h c A^  Is  presented  In  Conlon  (1977). 

In  this  transformed  state  space,  we  define  the  transition  probabilities 
to  be: 

P_(t)  • P[h(t)  » n|h(o)  - m] 

mn 

■ P[h(t  + s)  » n|h(8)  = m] 

■ V t + o(t), 

HUI 


?^(t)  - 1 - V t + o(t), 

m n 

where 


n*m 


and  h(t)  denotes  the  state  the  epidemic  Is  In  at  time  t.  The  hazard  rates 


for  the  epidemic  process  are 


V 

flBl 


n-kj^-k(8j^-l,  lj^+1;  Sj,  12) 
n“k2*k(Sj^,  lj^”l;  82'  ^2^ 
®2®2^1*  Sj-l,  12'*'^^ 

Y2l2»  n-k^-k(Sj^,  1^^;  s^.  12*1) 


A Markov  chain  Is  said  to  be  unlformlzable  If  sup  v ■■  v < «*.  Since 

m IB  O 
m 

our  chain  Is  finite,  unlforulzablllty  follows  Immediately.  Kellson  (1974) 
showed  that  the  original  process  (our  unlformlzable  chain)  and  a dlscrete-tlne 
process  constructed  by  randomly  selecting  the  nuadier  of  transitions  according 
to  a Poisson  process  with  parameter  have  state  probabilities  Identical 


in  law 
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For  soae  v 2 v , lot 
o 


V /v,  m * n 

■n 

1 - V /v,  m - n. 


ond 


K • <•„>■ 


Thon,  follovlng  Kollson  (1974),  we  can  derive  the  consistency  relation 


i»  - (T  ) - ilp-A,. 

■ V V VgO 

where  ^ Is  the  asitrlx  restricted  to  the  good  states. 


(3.1) 


To  Illustrate  the  theory,  let  us  consider  the  first  nine  states  of  an  epl- 
dealc,  as  presented  above.  If  we  choose  0j^*O.l,  Yj^~0.2,  and  Y2”0*2> 

then  the  are  given  by 

123456789 


0 0.2 


0.2  0.4  0.4 


0.3 


Thus,  ■ 


* 0.4.  Let  V ■ 1.  Then  the  matrix  A 


v-1 


Is 


1 0 

.2  .8 


0 

.2 

0 

0 

0 

0 

0 


0 

0 

.2 

0 

0 

0 

0 


0 

0 

1 

0 

0 

.2 

0 

0 

0 


0 

0 

0 

.8 


0 

0 

0 

0 


0 
0 
0 
0 

.2  .6  0 

0 .2  .6 

0 0 

.1 

0 


0 

0 


0 

0 

0 


0 

0 

0 

0 

0 

0 

1 


0 

0 

0 

0 

0 

0 

0 


.2  .7 

0 0 


0 

0 

0 

0 

0 

0 

0 

0 

1 


Ms  observe  that  G consists  of  the  state  corresponding  to  IifS  only.  Thus,  A^ 

a scalar,  and  the  consistency  equation  (3.1)  gives  the  naan  duration  tlna  of  an 
epldeale  which  Is  currently  In  state  (1,0;0,1): 


T,  . tt-0.7)-l.l  . 3.3. 


We  note  that  the  choice  of  v i Is  arbitrary.  If,  for  exaiiq>le,  we  choose 
V ■ 2,  we  have 

Tg  - i(l-.85)"^-l  - a.I. 

In  geaeral,  when  more  good  states  are  considered,  the  dimension  of  A _ 

V,G 

Is  quite  large,  and  the  Inversion  of  [Ig  - g]  Is  cumbersome.  For  example. 

In  the  situation  where  ■ N2^  ■ 8,  there  are  2025  states  In  the  transformed 

state  space  A^,  with  505  states  In  B and  1520  states  In  6.  A computationally 

feasible  algorithm  Is  required. 

4.  A RECURSION  APPROACH 

In  the  previous  section,  we  presented  the  consistency  equation  for  T , 

in 

together  with  the  relevant  theoretical  justifications.  We  also  noted  potential 

difficulties  In  actually  computing  values  for  T when  the  population  sizes  are 

m 

not  small.  We  conment  here  on  a computationally  feasible  (and  easily  program- 
mable) approach  to  the  calculation  of  T . 

n 

We  observe  that  the  mean  duration  time  of  our  epidemic  process  at  any 
given  state  can  be  decomposed  Into  the  sum  of  the  following  components:  the 
average  time  spent  In  the  given  state,  and  the  weighted  mean  duration  times 
of  the  epidemic  In  those  states  to  which  the  epidemic  may  proceed  In  one  step 
from  the  given  state.  The  weights  are  the  hazards  of  each  permissible  one- 
step  transition  from  the  given  state.  Thus,  If  we  let  T(s^,l^;s2>l2)  denote 


th«  mean  duration  time  of  an  epidemic  which  is  currently  In  state  (Sj^,lj^;s2tl2) 


we  may  write  the  following  recursion  relation: 


TCsjf  ij^5S2»i2) 


''^Sl*ll‘S2*V  f 


+ T(8j^,lj^-l;s2»i2) 

+ 82^2^1  T(82,ij;s2-l.l2+l) 


+ >2^2  T(8j^,l2;82.l2 


-1)  , 


where  v(s^,l2;82>l2)  ■ ^1*1^2  ^1^1  ^2®2^1  ^2^2' 


k computer  program,  based  on  this  recursion  relation,  has  been  written  to 
calculate  the  mean  duration  time  at  any  state  In  the  epidemic  process.  It  uses 
the  counting  system  described  In  Section  3 and  has  the  desirable  feature  that 
the  calculation  of  the  mean  duration  time  of  an  epidemic  In  state  h requires  the 
computation  of  mean  duration  times  for  only  those  states  j,  where  j < h.  CA 
Fortran  listing  of  the  program  Is  available  from  the  first  author  on  request.] 


5.  SCMB  NUMERICAL  RESULTS 

In  this  section  we  present  the  results  of  an  investigation  of  the 
relation  between  the  mean  duration  time  and  the  quantities  It  depends  on,  namely, 
the  infection  and  removal  rates  and  the  numbers  of  susceptlbles  and  Infectlves 
In  the  two  populations.  Our  investigation  restricts  attention  to  epidemics 
starting  In  a state  of  the  form  (Sj^,l:s2>l)>  l.e.,  with  one  initial  Infective  In 


each  population. 


- 10  - 


We  first  note  that  the  mean  duration  time  for  an  epidemic  starting  In 
state  (Sj^,l^;s2,l2)  with  parameter  vector  ^ times  larger 

than  the  mean  duration  time  for  an  epidemic  starting  In  the  same  state  with 
parameter  vector  (c3j^,cB2>‘cYj^>cy2^*  c > 0.  Hence,  only  one  parameter 
vector  in  the  family  { (cBj^,cB2;cYj^»CY2) » c ^ 0}  need  be  examined. 

For  parameter  vectors  of  the  form  (6, 6; 9, 6),  the  mean  duration  time  In- 
creases from  Its  initial  value  at  s^  ~ S2  * 1 to  a maximum  at  s^  ~ S2  * 3 of 

1.52  times  the  Initial  value  and  then  decreases  to  a value  at  s^  ■■  S2  ~ 8 

which  Is  1.05  times  the  Initial  value,  as  Illustrated  In  Figure  1. 

[Figure  1 about  here.] 

To  determine  why  the  mean  duration  time  Increases  and  then  decreases 
as  Sj^  * 82  increases,  the  parameter  vector  (0.5, 0.0; 0.0, 0.0)  was  investi- 
gated and  the  results  are  presented  In  Figure  2. 


[Figure  2 about  here.] 


r 


Figure  1 


1 
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Mean  duration  tine  for  epidemics  starting  In  state  (Sj^*  1;  S2*  1)  with 
parameter  vector  (0.5,  0.5;  0.5,  0.5). 


1 


MEAN  DURATION  TIME 
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Figure  2.  Mean  duration  tlae  for  epidemics  starting  In  state  (s^,  1;  S2>  1)  with 

i 

I parameter  vector  (0.5,  0.0;  0.0,  0.0). 
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For  this  parameter  configuration,  the  epidemic  can  end  only  if  Sj^  ■ 0, 

since  any  existing  Infectlves  are  not  removed.  The  consequence  of  this  is 
that  the  mean  duration  time  will  Increase  as  Sj^  ■ S2  increases  since  the 

state  (k,l;k,l)  decays  into  (k-l,2;k,l),  which  has  Che  same  mean  duration 
time  as  (k-l,l;k-l,l)  for  these  parameters. 

Examination  of  Figure  3,  which  presents  Che  mean  duration  Claw  for  an 
epidemic  with  parameter  vector  (0.5,0.5;0.0,0.0),  indicates  that  allowing 
susceptible  carriers  to  become  Infected  (B^  * 0)  results  in  a decrease  in 
the  mean  duration  time,  for  sufficiently  large  s^  • *2*  the  con- 

tinued Increase  as  in  Figure  2. 

[Figure  3 about  here.] 

With  this  parameter  configuration,  the  epidemic  may  terminate  only  if 

s^  * 0.  The  decrease  in  mean  duration  time  is  due  to  the  fact  that  the  rate  of 

ttansltlons  of  type  (1)  is  B2S2I2*  Since  B2  * 0,  I2  Increases  and  the  transition 

rate  tends  to  be  larger  than  it  is  in  the  situation  in  Figure  2,  where  I2  remains 

constant.  Allowing  nonzero  removal  rates  does  not  affect  the  general  shape  of 
the  mean  duration  time  curve,  as  can  be  seen  by  comparing  Figures  1 and  3. 

Thus,  the  Increase  and  subsequent  decrease  in  the  isean  duration  time  ob- 
served in  Figure  1 can  be  explained  as  follows.  An  epidemic  starting  in  state 
(1,1;1,1)  with  parameter  vector  (6, 6:6, 6)  is  most  likely  to  terminate  because  s^ 

becomes  zero.  As  another  host  and  carrier  are  added,  the  time  required  for  Sj^ 

to  decrease  to  zero  Increases  and  the  epidesilc  lasts  longer.  However,  the 
addition  of  the  fourth  individuals  causes  an  Increase  in  the  hazard  rates. 


MEAN  DURATION  TIME 


which  More  than  offsets  the  additional  time  required  to  remove  the  fourth 
susceptible  host.  Thus,  the  mean  duration  time  of  an  epidemic  In  state 
(4, 1:4,1)  la  lass  than  that  of  an  epidemic  In  state  (3,1; 3,1). 
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20.  ABSTRACT 

In  this  paper,  the  two-population  model  for  a carrier-borne  epidemic  posed  by  Bailey 
(The  Mathematical  TOeory  of  Infectious  Diseases  and  Its  Applications.  1975,  p.  211)  is 
foraulatad  In  a mathematically  tractable  manner.  This  sndel  reflects  the  epidemiology 
of  diseases  such  as  malaria,  where  the  progress  of  the  disease  depends  on  the  Interaction 
of  a population  of  mosquitoes  and  a population  of  humans.  An  expression  for  the  mean 
duration  tlM  of  the  epidemic  Is  obtained  and  a computationally  feasible  algorithm  is 
presented.  Results  of  a study  Investigating  the  consequences  on  the  mean  duration  time 
of  varying  the  Infection  and  removal  rates  In  the  two  populations  are  given. 


